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We give a theoretical study of a double-cavity system in which a mechanical resonator beam is coupled to 
two cavity modes on both sides through radiation pressures. The indirect coupling between the cavities via the 
resonator sets up a correlation in the optomechanical entanglements between the two cavities with the common 
resonator. This correlation initiates an entanglement transfer from the intracavity photon-phonon entanglements 
to an intercavity photon-photon entanglement. Using numerical solutions, we show two distinct regimes of 
the optomechanical system, in which the indirect entanglement either builds up and eventually saturates or 
undergoes a death-and-revival cycle, after a time lapse for initiating the cooperative motion of the left and right 
cavity modes. 


I. INTRODUCTION 

Cavity optomechanical systems (Hi arise from the classi¬ 
cal Fabry-Perot interferometer [j2[] by replacing one of the 
fixed sidewalls with a cantilever or double-clamped beam [3|- 
U]. The one-dimensional degree of freedom introduced by 
the movable mechanical element adds a free resonator mode 
to the cavity system and allows this mode to interact with 
the cavity field through radiation pressure on the reflectively 
coated mechanical resonator. Regarded as a micromirror, this 
resonator can be feedback-controlled through the cavity field, 
on which numerous cooling protocols have been conceived 
and experimentally demonstrated in the last decade IM3- 

The degree of control in this hybrid cavity-micromirror sys¬ 
tem can be further enhanced when the micromirror is replaced 
by a double-face reflective membrane jT 1 , fl2ll . If a second op¬ 
tical cavity is coupled to it on the opposite side of the existing 
cavity, a two-mode or double-cavity optomechanical system 
with enhanced nonlinearity is formed ll3lTl6ll . Entanglement- 
wise, though it was observed that the enhanced squeezing re¬ 
sulted from the nonlinear coupling helps generate static entan¬ 
gled state of distant mirrors [130, the dynamic property of en¬ 
tanglement between the two cavities is less well-understood. 

Recent studies reveal that the dynamics of phonon-photon 
entanglement plays an important role in defining the sys¬ 
tem characteristics, such as the transitions between oscillation 
modes mm, robustness against noisy environment m, 
sudden death and revival of states I120l422ll . and optimal entan¬ 
glement |[23f l. In this article, we study the dynamics of the en¬ 
tanglements in a double-cavity optomechanical system where 
each photon mode in the two opposite cavities is, structure- 
wise, symmetrically coupled to a common mechanical res¬ 
onator mode via radiation pressures, albeit asymmetric cou¬ 
pling strengths and driving powers are generally assumed. 
Our main concern is to determine how the cavity-resonator 
entanglements [24] can be transferred to the indirectly cou¬ 
pled cavities over time. 

We show here such an entanglement transfer is possible in a 
double-cavity optomechanical system through measuring the 
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entanglements in logarithmic negativities among the compo¬ 
nent pairs. In particular, the negativity is computed through 
determining the symplectic eigenvalues of a covariance ma¬ 
trix that relates the fluctuations of all six quadratures of the 
system’s main components. This method is standard in the 
literature of dynamic entanglement but we have generalized it 
to apply on a 6 x 6 covariance matrix. We observe that the 
successful generation of entanglement transfer only requires 
a single-sided driving laser and that the transfer patterns can 
be distinctively categorized into two groups for the different 
operating regimes assumed by optomechanical system. 

Moreover, all the logarithmic negativities computed exhibit 
a time delay before the first appearance of a non-zero value. 
This time point signifies the initiation of cooperative motions 
among the three components in the optomechanical system, 
showing the transient response of the system to the external 
driving lasers as a whole. Nonetheless, the indirect entangle¬ 
ment between the left and the right cavity is apparent in all 
cases, thereby facilitating a mechanism for entanglement re¬ 
lay through cascaded cavities although the cavities are physi¬ 
cally not directly coupled. Such a mechanism would be use¬ 
ful to quantum information processin g, espe cially in terms of 
non-adiabatic quantum state transfer 112 5l 12611 . and would pro¬ 
vide a physical means to realize cavity arrays or resonator 
waveguides for transmitting information encoded in a quan¬ 
tum state l27l - l29ll . 

In Sec. QT| we give a detailed description of the double¬ 
cavity model. The equations of motions are derived under the 
Heisenberg picture in Sec. [Ill] and the steady-state solutions 
are calculated to give proof of the sufficiency of single-sided 
driving. After the covariance matrix of the fluctuations is in¬ 
troduced, the entanglements among all component pairs are 
computed numerically and analyzed in Sec. [IV] The conclu¬ 
sions are given finally in Sec.lVl 


II. DOUBLE OPTOMECHANICAL CAVITY 

The proposed double-cavity optomechanical system is il¬ 
lustrated in Fig. Q] in which a mechanical resonator with re¬ 
flective coatings on both sides receives the radiation pressures 
from both the cavity on the left side (L) and the cavity on the 
right side (R). The total Hamiltonian H = Hq + H rAr \ + H exi 
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FIG. 1. (Color online) Model schematic of the double-cavity op¬ 
tomechanical system: a mechanical element with reflective coatings 
on both sides serves as a double-face mirror that experiences radia¬ 
tion pressures from both the left-side cavity and the right-side cavity. 
An incident driving laser enters the double-cavity system from the 
left side. 

thus consists of three parts, which reads (h = 1), respectively, 

1 

H 0 = u L a' L a L + LO R a^ R a R + — + -mVL 2 M q 2 , (1) 

#rad = ( 'vLa^dL - g R a^ R a R ^ q , (2) 

#ext = i^L {a\e~ luld ' Lt - h.c.) + ie R 

x (a^e~ iuid ’ Rt - h.c.) (3) 

The part H 0 accounts for the free Hamiltonians of the res¬ 
onator and the cavities, the latter being regarded as bosonic 
modes of frequencies u>l and ui R . We associate a pair of an¬ 
nihilation and creation operators a a and ag f or eac h bosonic 
mode, where a indexes the cavity side, either left L or right 
R. We assume the frequencies ujl and cj r to be different in 
general according to the asymmetric cavity lengths t a and fi¬ 
nesses F a assumed. The leakage rates are defined correspond¬ 
ingly from these parameters: n„ = irc/2F a £<j. The mechani¬ 
cal resonator is described by the conjugate pair q and p, along 
with its oscillation mode frequency of and its mechanical 
damping rate of Tm- 

The part // ra( | accounts for the phonon-photon interactions 
derived from the radiation pressures. The radiation pressure 
from each cavity side results from the deformation of the cav¬ 
ity volume due to the displacement q of the middle resonator, 
which shifts the resonance frequencies of each cavity modes. 
Expending this frequency shift to first order, i.e. g /, ~ wl/ 4, 
the radiation pressure term sets the optomechanical coupling 
on the left side with strength p L = g L /yJ2vnF\M, where m 
is the effective mass of the resonator mode. The same deriva¬ 
tion applies to the right cavity, giving the coupling strength 
Vr = 9r/ y/2mQ,M but the leading sign would be opposite 
to ///, as the common resonator has its displacement q follow 
opposite directions for the two radiation pressures. Between 
the two cavities, there is no direct coupling. 

The part H ext accounts for the two external driving lasers 
with frequency uid,L and frequency uJd,R■ The driving strength 
s& of each laser is related to the input laser power P a and the 
leakage by \e a \ 2 = 2hi a P a /huid,a- Note that we assume an 
asymmetric setting for the double-cavity system: the radiation 
pressures from the two sides are not identical and the cavities 
are unequally driven. 


To study the indirect entanglement across the two cavity 
modes, we begin with the dynamics of the three components 
in the double-cavity optomechanical system through deriv¬ 
ing a set of nonlinear Langevin equations. We carry out this 
step by finding the Heisenberg equations of the operators from 
the Hamiltonian in E. ([Hi-® and introducing phenomenologi¬ 
cally the relaxation terms and their associative Brownian noise 
terms. The Langevin equations under the rotating frames of 
reference read 



p = - mP? M q - T mV- riL(J L a L + r/ R a^ R a R + 

CiL = — («L + iA L )a L - ir]Ld L q + £l + \/2«lo“, 
clr = ~(kr + iA R )a R + ip R a R q + e R + \Z2k r o, r , (4) 

where A 0 ,l = oJl ~ u d ,L (A 0 ,fi = u R - u d ,R ) is the static 
detuning of the left (right) cavity field from the left (right) 
driving laser. The zero-mean fluctuation terms obey the 
correlation relation (a“(f)a^(f')) = S(t — t 1 ). 

The mechanical mode is under the influence of stochastic 
Brownian noise that satisfies in general the non-Markovian 
auto-correlation relation with a colored spectrum: 

S {»* h ) + 

(5) 

where k R is the Boltzmann constant and T is the temperature 
of the mechanical bath. However, for a high quality mechan¬ 
ical resonator with Q = G .\,f /V M 1, this non-Markovian 
process can be approximated as a Markovian one, where its 
fluctuation-dissipation relation can be asymptotically simpli¬ 
fied to 0 [H : 

man+atm) /2 = r M (2n +mt - o, (6) 

where n = (exp{ VLm/ ksT} — l) -1 is the mean occupation 
number of the mechanical mode. This simplified Markovian 
relation will be assumed in the calculation of the entangle¬ 
ments. 


HI. DYNAMICS AND ENTANGLEMENT 
A. Steady states 

In a single-cavity optomechanical system, the radiation 
pressure contributes the nonlinear photon number term in the 
Langevin equation of the mirror momentum p in Eq. (0), lead¬ 
ing to a multistability of the coordinate p with three nonzero 
steady states. For a double-cavity case here, the second radia¬ 
tion pressure by the other cavity contributes a similar term in 
the equation. Under the asymmetric setting, the two pressure 
terms are not commensurate and the number of steady states 
of p increases to five. The steady states are given by 

/ a _ ~Vl \{a L )\ 2 + m |(aii )| 2 ^ 

W mfljj ’ U) 

^ k ct +*(A o,aAr] a {q)y ^ 
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where the plus (minus) signs in the second equation refers to 
the left (right) cavity. 

For entanglement generation, it is necessary for the equa¬ 
tion set 0-0 to have non-zero steady states. Therefore, a 
single-cavity optomechanical system usually requires an ex¬ 
ternal driving laser (i.e., non-zero value of e) to drive the me¬ 
chanical resonator out of its zero steady states at equilibrium 
position. However, for optomechanical systems with double¬ 
sided cavities, one external driving laser at either end of the 
cavities is sufficient to drive the mechanical resonator out of 
its zero position, in which case Eq. 0 would fall back to the 
single-cavity case of three roots. 

In addition, we can observe that even when the double cav¬ 
ities have exactly symmetrical setup, i.e. identical laser driv¬ 
ing amplitudes (£l = er = e), radiation pressures ( 77 ^ = 
VR = V)’ and cavity lengths, the differing signs before r] a ( q) 
to be taken by (cll) and (cm) in Eq. 0 allows the cavities 
to admit non-zero steady states. This is because the two cav¬ 
ity modes are constructively interfering with each other at the 
interface of the mechanical resonator through their indirect in¬ 
teractions of radiation pressures. In other words, even though 
the radiation pressures are exerted along opposite directions, 
the dynamic 7r-phase difference between the cavities fields, re¬ 
flected in the Hamiltonian Eq. 0 as the generator of the cav¬ 
ity motion, render the radiation pressures out of phase to favor 
the generation of entanglement. Given the symmetric setting 
where kl = Rr = n and Ao,l = Ao,k = Ao in addition 
to the identities in driving amplitudes and radiation pressures, 
the condition for the steady-state equations to admit real roots 
is the inequality among the system parameters 


, ' £i Vw ( '‘ 2+A » ) - (9) 

Its derivation is given in Appendix A. Finding m£l 2 M as the 
Young’s modulus of the resonator < Ao) and that the 

cavities have sufficient finesses (n < Ao), the above crite¬ 
rion is met in most scenarios and the validity of entanglement 
generation is almost guaranteed. 

For the symmetric setting, we expect the patterns of entan¬ 
glement generations between either end of the cavity modes 
and the mechanical resonator to be qualitatively similar and 
differ only quantitatively in their variations over time. Devi¬ 
ating from this setting, the increase in asymmetry among the 
system parameters would increase the qualitative difference 
between the patterns of entanglements. We demonstrate these 
effect later in Sec.HVl 


B. Entanglement measure 


Theoretically, the entanglements in terms of logarithmic 
negativity are computed through the fluctuations of the cav¬ 
ity quadratures about the steady states obtained from Eqs. 0- 
®. That is, we define the dimensionless quadratures of the 


two cavity fields as 


Y a 


71 (° CT + °t) ’ 


GO) 

(ii) 


and the corresponding input noise operators accordingly. 
Then taking O = (q. p. X R , Y[ . X R , Y r ) as the vector oper¬ 
ator for all the quadratures in the optomechanical system, we 
expand it to first-order using a c-number steady-state value 
and a zero-mean fluctuation operator 0(t) = ( O) + 60(t). In 
addition, the nonlinear terms are linearized assuming | (a) | 

1 in the expansion: (ala) ~ | (a) | 2 and (aq) — (a) ( q ), while 
the higher-order products of the fluctuations are ignored. 

The Langevin equations in Eq. 0 with the first-order ex¬ 
pansion gives a coupled system of differential equations about 
the noise operators, enabling the coupling between the fluc¬ 
tuations of the two cavity fields and the mechanical resonator 
and thus the generation of entanglement between the two opti¬ 
cal modes. Note that even though we have linearized the equa¬ 
tions for these operators, eliminating the mechanical quadra¬ 
tures q and p in Eq. © will lead to equations of and a' L 
nonlinearly related to a R and a} R . This implies that the indi¬ 
rect entanglement between the quadratures of the left and the 
right cavities follow a nonlinear form in time. 

In the following, instead of solving the coupled equa¬ 
tions analytically, we follow the standard numerical ap¬ 
proach adopted by the current researches on dynamic entan¬ 
glement (17, [l_8l [24J]. The difference here is that we have 
a 6-component vector u = (Sq, Sp. dXr,, SYj,, SX R . 6Y R ) 
over the six quadratures of the tripartite optomechanical sys¬ 
tem instead of the usual 4-component vector. Similarly 
extending the input-noise vector to the 6-component n = 
(0, £, y/2 k l X£, sJ2k l YE 1 , y/2K R X'g, yJ2n R Y£), we write 
the time-dependent inhomogeneous equations of motion as 
u{t) = A(t)u(t) + n(t), where A(t) = 


0 1/m 0 

-mVL 2 M -r M -G xL (t) 

Gy L (t) 0 -k l 

—G X L,(t ) 0 —A L (t) 
-Gy R {t) 0 0 

. G xR {t) 0 0 


0 0 0 " 

— Gyl,{t) G xR (t) Gy R (t) 

A L (f) 0 0 

—kl 0 0 

0 — hr A R (t) 

0 —A R (t) — kr _ 


( 12 ) 

In the matrix, G xa {t) = r] a (x(t)} and G y<J (t) = r) a (y(t)} 

are the real and the imaginary parts of the scaled coupling 
constants G a (t) = V2(a a (t))r) cr . Along with the oscillation 
of the mechanical resonator, the dynamic detunings of the two 
cavities are defined as 


A ff (f) = Ao.^ ±r] a {q(t)), (13) 

where the plus (minus) sign corresponds to the left (right) cav¬ 
ity. 

When the tripartite system is stable, it reaches a unique 
steady state, independently from the initial condition. Then 
given any arbitrary steady state, the fluctuations about it is 
fully characterized by its 6 x 6 covariance matrix V of the 
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pairwise correlations among the quadratures, which obeys the 
equation V(t) = A(t)V(t)+V(t)A T (t)+D. The diagonal el¬ 
ements of the V are, in order, auto-correlations of the quadra¬ 
tures of the resonator, the left, and the right cavity mode. 
Hence, D = diag(0, Tm( 2n+l), kl, kl, kr, Rr) is the diag¬ 
onal matrix for the corresponding damping and leakage rates 
responsible for the fluctuations. The multiple fluctuation- 
dissipation relations defined in Sec.Q]]are therefore encapsu¬ 
lated in the relation (rii(t)nj(t') + n,j(t')ni(t )} /2 = S(t — 
t')Dij. From its evolution equation, the covariance matrix V 
can be written as a block-matrix 


V = 


Vm Vml 
Vml Vl 
[ V mr Vlr 


Vmr 

Vlr 

Vr 


(14) 


where each block represents 2x2 matrix. The blocks on the 
diagonal indicate the variance within each subsystem (the res¬ 
onator M, the left cavity mode L, and the right cavity mode 
R), while the off-diagonal blocks indicate covariance across 
different subsystems, i.e. the correlations between two com¬ 
ponents that describe their entanglement property. 

To compute the pairwise entanglements, we reduce the 6 x 6 
covariance matrix V to a 4 x 4 submatrix Vs. There are three 
such cases of the submatrix Vs: (i) if the indices i and j for the 
element V. l:l are confined to the set {1, 2, 3,4}, the submatrix 
Vs = [Vij] is formed by the first four rows and columns of V 
and corresponds to the covariance between the resonator mode 
and the left cavity mode. Similarly, (ii) if the indices run over 
{1, 2, 5, 6 }, Vs is the covariance matrix of the resonator and 
the right cavity mode, (iii) If the indices run over {3,4, 5, 6 }, 
Vs designates the covariance between the two opposite cavity 
modes. Summarizing, the submatrix can be written as 


Vs 


V Q 

vL 


Vx /3 

Vp 


(15) 


where a, /?, and 7 index the subsystems {M, L, R} in the op¬ 
tomechanical cavity. The entanglement measured by logarith¬ 
mic negativity is computed through a process known as sym- 
plectic diagonalization of each submatrix Vs, where the en¬ 
tanglement properties are contained in the symplectic eigen¬ 
values of the diagonalized matrix. If we write the diagonalized 
matrix as diagfu , v , u_i_, then the eigenvalues along 
the diagonal read [32] 


1 r 


77 S(Vs)T ^(Vs ) 2 — 4 detVs 


2 L 


(16) 


where £(Vs) = det(V„) + det(Vg) — 2 det(V 1| g). 

Denoting the state of a bipartite subsystem in the tripartite 
optomechanical cavity as p , the negativity is defined as 

||p T |L -1 

N(p) = -, (17) 

where || p T || 1 indicates the trace norm of the partial transpo¬ 
sition of p 13311 . Taking u_ as the minimum symplectic eigen¬ 
value of the covariance matrix, ||/9 T || 1 is equivalent to l/u_ 



FIG. 2. (Color online) Time evolution of the tripartite optomechani¬ 
cal system characterized by the logarithmic negativities between (a) 
the left cavity mode and the mechanical resonator, (b) the right cav¬ 
ity mode and the mechanical resonator, and (c) the left and the right 
cavity modes. Two cases are shown with different colors: (i) the blue 
curves show the symmetric case where the parameters of the left and 
the right cavities are set identical; and, in contrast, (ii) the red curves 
show the asymmetric case where some parameters of the two cavities 
are set distinct. The parameter values taken for the plots are given in 
the text. 


after the diagonalization. Hence, the negativity is a decreas¬ 
ing function of and we usually write N(p) = max{0, (1 — 
) / 2 u_ } and take its logarithmic value Em = In 11 p T 11 as a 
measure of the entanglement [34ll. This logarithmic negativity 
has the expression Ejy = max{ 0 , — ln(u_)}. 

In other words, the symplectic eigenvalue V- completely 
quantifies the quantum entanglement between each pair of 
components in the system. The necessary condition for show¬ 
ing a bipartite subsystem is entangled is that the symplectic 
eigenvalue retains a value less than one, which is equivalent 
to the inequality 4det Vs < E(Vs) — 1/4 ll35ll . 


IV. ENTANGLEMENT TRANSFER 
A. Delayed build-up 

To measure the entanglements, the noise terms £, A/ 1 , Y[ n , 
X , and V / 1 ) 1 that appear in the variance matrix of Eq. (IT4l) 
are taken as random variables of zero-mean Gaussian pro¬ 
cesses. The entanglements measured in logarithmic negativi¬ 
ties Em are plotted against time for each of submatrices given 
in Eq. i | 13 [ 1 to discern the entanglement transfer. We found 
similar transfer patterns over a range of parameters close to 
the experiments 1361. One typical case is shown here in Fig. [2] 
where from top to bottom we plot, respectively, Em between 
the left cavity and the resonator, between the right cavity and 
the resonator, and finally between the left and the right cavi¬ 
ties. 
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For comparison, two cases are plotted for each entangle¬ 
ment pair: the blue ones denote the symmetric case and the 
red ones denote the asymmetric case. For the symmetric 
case, we adopt for the mechanical resonator a quality factor 
Q = 20000, resonance frequency u>m = 1MHZ, and effective 
mass m = 10 ng ; for the cavities, we take cavity length at 22 
mm with finesse F = 2.6 x 10 5 and cavity mode wavelength 
of 1064 nm. We set the power of the driving lasers at 70//W, 
which is detuned from the cavity mode at A = 6.5 u>m- For 
the asymmetric case plotted in red, we have adjusted the right 
cavity to a length of 19 mm, which consequently affects the 
cavity leakage and the coupling amplitude between the driv¬ 
ing and the cavity, while the length of the left cavity and other 
parameters remain unchanged. 

We observe from Fig. [2]that there are two phases in the en¬ 
tanglement evolution. The initial phase is a period of zero E jv, 
showing a delay in the formation of entanglement. The latter 
phase is a gradual build-up until certain saturation is reached. 
While the entanglement generations between either cavity and 
the mechanical resonator are smooth, that between the two 
cavities are oscillating or quasi-oscillating because of the non¬ 
linear nature of the radiation pressure coupling 01. Averag¬ 
ing out the oscillation, we see the patterns in the build-up of 
entanglement are identical to those between the cavity and the 
resonator. In addition, the delay periods among all three pairs 
coincide, demonstrating the transfer of cavity-resonator entan¬ 
glement to intercavity entanglement and showing that distant 
entanglement is possible if the distant objects are indirectly 
coupled. 

The delay in the entanglement build-up, during which Ejr 
assumes zero value, corresponds to the negativity in Eq. 
taking a nonphysical negative value. We can interpret this de¬ 
lay period as the time duration when the three components 
in the tripartite system spend to establish their cooperation, 
which like the effect of superradiance depends strongly on the 
resonance linewidths. Comparing the delays for the symmet¬ 
ric and the asymmetric cases from Fig.[2 a ) and (b), we see 
the similar inverse proportionality in the entanglement delay 
Tr> on the cavity leakage rates K a , i.e., Tr> oc ft” 1 . When 
the cavities are setup symmetrically, we measure the delays 
in both Fig. [2 a) and (b) at about 89/xs; when they are setup 
asymmetrically with hr < hr, we observe Tr> for the left 
cavity being greater than its counterpart at the right side, at a 
difference of 15.7/iS in time for a difference about 2.3kHz in 
cavity linewidths. 


B. Death and revival 

The influences of asymmetric parameter setup for the cav¬ 
ities are not only reflected in the delays of entanglement gen¬ 
eration, but also in the entanglement pattern itself. In Fig. [3j 
we show a typical example with entanglements generated in a 
pattern distinctly differently from those in Fig. [2] The entan¬ 
glements measured in logarithmic negativity are again plot¬ 
ted from top to bottom, respectively, for the three component 
pairs discussed above, but with driving laser powers increased 
to 80/uW and cavity finesses decreased to F = 1.0 x 10 5 . The 
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0 0.4 0.8 1.2 1.6 2.0 


Time (xl0“ 4 s) 

FIG. 3. Time evolution of the logarithmic negativity En for the same 
three pairs of components in the tripartite system as in Fig.[2]plotted 
respectively in (a), (b), and (c), where (c) shows the death and revival 
patterns in the intercavity entanglement. 


left and the right cavity lengths remain in an asymmetric setup 
of 22mm and 20mm, respectively, and the rest of parameters 
are kept identical to those in Sec. II VI A 

While the cavity-resonator entanglements for the two cav¬ 
ities follow the pattern of build-up to saturation after a time 
delay, which is similar to those of Fig. [2 the intercavity entan¬ 
glement does not but otherwise oscillate over a death-revival 
cycle. Because of the inverse proportionality of the time de¬ 
lay to the cavity linewidths, the plots show a shortened delay 
and a reduced discrepancy between the delays in the left and 
the right cavity-resonator entanglements due to the decrease 
in cavity finesses. 

On closer inspection, we can see the build-up in (a) and 
(b) are sharper and less gradual than their counterparts in 
Sec. ll VT \ and the absolute negativity they can obtain are much 
smaller, especially for the left cavity. Even for the right cav¬ 
ity, its entanglement with the resonator declines shortly after a 
peak value, making all three plots assume essentially different 
characteristics than those of Fig. [2 This distinction can be at¬ 
tributed to the strong dependence of the operating regimes of 
optomechanical systems on external driving power and cavity 
finesse. In a single optomechanical cavity, it is reflected as 
periodic and quasiperiodic motions of the resonator JH ; in 
the double optomechanical cavity here, it is reflected as the 
resonator being driven monotonically in-phase (Fig.|2 c )) and 
driven periodically in-phase and out-of-phase (Fig. [2c)) with 
the left and right cavities. 


V. CONCLUSIONS 

To summarize, we have studied the dynamic transfer 
of quantum entanglement from those within two cavity- 
resonator pairs to that between these two cavities inside a 
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double-cavity optomechanical system. We numerically solved 
a coupled set of Heisenberg-Langevin equations to show the 
generation of quantum entanglements between each pair of 
the components under an experimentally accessible set of pa¬ 
rameters. We find that the entanglement of the indirectly cou¬ 
pled cavities is built up over time in a pattern similar to those 
of the directly entangled cavity-resonator pairs, verifying the 
entanglement transfer. The similarities are accentuated by the 
almost identical characteristic delays and rising patterns but 
the entanglement transfer would be suppressed by the asym¬ 
metries in the two cavities. The asymmetries also differenti¬ 
ates the initiation times of the cavity-resonator entanglements, 
which leads to our speculation that the tripartite system is un¬ 
dergoing a cooperation process similar to that of superradi¬ 
ance before the emergence of the entanglement. To under¬ 
stand such a transient effect in a multipartite system demands 
a detailed analysis of the Heisenberg-Langevin equation set, 
which we shall leave to future studies, but we have seen here 
that dynamic entanglement is not only a measure of quantum 
information, but also a useful tool to dissect the cooperative 
motions of microscopic systems. 
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Appendix A: Steady states of symmetrical double cavity 
optomechanical system 

Substituting Eq. © into Eq. CZ} and cancelling the factor 
(q) on both sides of the equation, which implies the trivial 


solution being one of the steady state in the symmetrical cavity 
setup, we arrive at the quartic equation 





4Aq£ 1 2 3 4 5 6 7 

rmfn 2 M 


= 0 . 
(Al) 


Lacking the odd-order terms in ( q ), the roots (<y) 2 of the equa¬ 
tion can be solved directly through quadratic formula. Since 
n 2 + Ag > 0, the real roots (q) exist only when: 


i) (q) 2 is real, i.e. the discriminant being non-negative, 
which gives 


(rje) 2 > mfl 2 M K 2 A 0 , (A2) 

and ii) the quadratic root (q) 2 to Eq. (I A1 b is non-negative. 

To satisfy the latter, we have to consider two cases: 

ii-1) when k 2 — Aq < 0 , the square root of the determinant 
could take either the positive or the negative value. For the 
negative case, it is required that 

(ri e) 2 < (A3) 

or ii-2) for the postive case or when t; 2 — Aq > 0 , it is required 
that 

(qe) 2 > mn 2 M ^ . (A4) 

When two cases of condition (ii) are combined with con¬ 
dition (i), we see case (ii-1) impose a very stringent con¬ 
straint on the admissible values of {qe) 2 '. between zero and 
mtt 2 M (n 2 — Aq) 2 /4Ao- Case (ii-2) is more inclusive, which 
is what we are interested in here. Since it always holds that 
(k 2 + Aq) > 4k 2 Aq, when the inequality of Eq. (I A 41 holds, 
the first condition in Eq. (IA2b is automatically satisfied. 

To simplify the study, we confine our investigation in the 
positive domain of the detuning Ao, for which Eq. ( |A4 | i can 
be further reduced to 


qe > 



(« 2 + A 2 ). 


(A5) 
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